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ABSTRACT 



Aims. We analyze the carbon monoxide emission around the star WR 16 aiming to study the physical characteristics of the molecular 
gas linked to the star and to achieve a better understanding of the interaction between massive stars with their surroundings. 
Methods. We study the molecular gas in a region ~ 86' .4 x 86' .4 in size using CO J=l-^0 and '^^CO J=l— >0 line data obtained with 
the 4-m NANTEN telescope. Radio continuum archival data at 4.85 GHz, obtained from the Parkes-MIT-NRAO Southern Radio 
Survey, are also analyzed to account for the ionized gas. Available IRAS (HIRES) 60 yum and 100 yum images are used to study the 
characteristics of the dust around the star 

Results. Our new CO and '^CO data allow the low/intermediate density molecular gas surrounding the WR nebula to be completely 
mapped. We report two molecular features at -5 km s"' and -8.5 km s"' (component 1 and component 2, respectively) having a good 
morphological resemblance with the Ho- emission of the ring nebula. Component 2 seems to be associated with the external ring, 
whilst component 1 is placed at the interface between component 2 and the Ho- emission. We also report a third molecular feature 
~ 10' in size (component 3) at a velocity of ~9.5 km s"' having a good morphological correspondence with the inner optical and IR 
emission, although high resolution observations are recommended to confirm its existence. The stratified morphology and kinematics 
of the molecular gas could be associated to shock fronts and high mass-loss events related to different evolutive phases of the WR 
star, which have acted upon the surrounding circumstellar molecular gas. An analysis of the mass of component 1 suggests that this 
feature is composed by swept-up interstellar gas and is probably enriched by molecular ejecta. The direction of the proper motion of 
WR 16 suggests that the morphology observed at infrared, optical, radio continuum, and probably molecular emission of the inner 
ring nebula is induced by the stellar motion. 

Key words. ISM: molecules, radio continuum, ISM: Ring nebulae. Individual object: Stars: WR 16, HD 86161 
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1. Introduction 



The interstellar medium (ISM) surrounding massive stars is 
expected to be strongly modified and disturbed by the stel- 
lar winds and ionizing radiation. There is now consensus that 
massive stars (> 30 Mq) evolve trough a sequence of three 
stages, each of them characterize d by a different kind of wind 
dGarcia-Segura & Mac Lowlll995a b/). During the long-lived O 
phase (f ~ (2 - 4) X 10^ yr), the gas around the star is first ion- 
ized by the high Lyman continuum flux, producing an Hii region 
that is later evacuated by the powerful stellar w inds creating 
an "interstellar bub ble" (IB) (ICastor et ai]ll975t lOvson 1977,; 
I Weaver et al.lll977l ). The latter have been successfully detected 
both as thermal radio continuum shells and in the 21 cm line 
of at omic hydrogen (Arnal 1992; Cappa et al. 1996;l Arnal et alJ 
11999 ) in the form of shells having diameters of tens of parsecs, 
and expansion velocities of 6 - 20 km s ' . In the short-lived red 
supergiant (RSG) (f ~ (2 - 3) x 10^ yr) and luminous blue vari- 
able (LBV) {t ~ 1 X lO'* yr) stages, the stellar wind becomes 
denser and slower and ejections of up to tens of solar masses 
take place. The Wolf-Rayet (WR) stage (f ~ (2 - 6) x 10^ yr) is 
the last evolutionary phase of the massive stars prior to its ex- 
plosion as a supernova, and is characterized by a fast and chem- 
ically enriched stellar wind which rapidly reaches and inte racts 
with t he previous RSG/LB V wind creati ng a ring nebula (IChul 
[TgsTtlChu & Treff^ersll981allbl;IChu et al.i ri982). Hence, the mu- 
tual interactions between stellar winds and ejecta of different 



stages may form multiple shells in the ISM at scales from tenths 
to tens of parsecs during the star evolution dOarcia-Segura et sl\ 
fT996a"S). 

Since studies of the ISM around WR stars can provide im- 
portant information about the interaction between massive stars 
and their environs, molecular line observations are crucial instru- 
ments to understand that interaction. The molecular gas around 
WR stars has been analyzed in the last few years. One of the best 
st udied cases is the ri ng nebula NGC 2359 around the star WR 
7 (iRizzo et al.ll2001allbE ICappa et alJl200lh which shows clear 
signs of interaction between the star and its molecular and ion- 
ized environment. The stratification of the kinematics, morphol- 
ogy, and density of the molecular gas with respect to the central 
star can be associated to several energetic events related with dif- 
ferent evolutionary s tages previous to t he actual WR phase like 
RSG or LVB phases ('Rizzo et alJl2003 l). 

The Ha ring nebula around the WN8h star WR 16 (=HD 
86161) d van de r Hucht 2001), was first noticed bv lMarston et al.l 
(11994 . The nebula consists of a double ring system with a very 
intense inner ring ~ 5' north from the star, whilst the outer ring 
is placed to the northwest at ~8'. The emission of the outer ring 
appears to extend to the southeast of the star For the sake of 
clarity, in Fig[T]the Ho' emission distribution around WR 16 ob- 
tained from the Super Cosmos Ha Survey (SHSflis shown. The 
stratified distribution of the Ha emission around WR 16 may be 



http://www-wfau.roe.ac.uk/sss/halpha/index.html 
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Right Ascension (J2000) 



Fig. 1. Super Cosmos Hor image of the ring nebula around WR 
16 



indicative of an evolution of the ISM associated with different 
stages prior to the WR phase, which makes this nebula an excel- 
lent laboratory to study the interaction between WR stars with 
their environs. iMarston et alJ (119991) surveyed the COJ=l— >0 
line in a grid of ~20'x 14' centered at WR16 using SEST (HPBW 
= 45" ). They found a molecular "cocoon" at a radial velocitjQ 

of 5 km s"' around WR 16. They estimated a total molecular 

mass in the range 5-78 Mq. This fact, together with a chemical 
enrichment observed towards the inner part of the nebula, led the 
authors to suggest an stellar origin for both, the inner nebula and 
the molecular gas. 

Spectrophotometric distanc es of WR 16 were e stimated 
by several authors: 2.6 kpc (fConti & Vacca' '199(1, 2 kpc 
(IMarston et al.ll 19991) . 2.37 kpc (van der Hucht 2001). In this pa- 
per we adopt an averaged value of 2.3 ± 0.3 kpc. 

Although providing important information about the molec- 
ular gas associated with the close surroundings around WR 16, 
the observations of Marsto n et al.l (I1999D do not provide an en- 
tire picture of the molecular environment of the nebula because 
they observed the distribution of the molecular gas in a limited 
region around the WR star. Clearly, the CO emission distribution 
extends well beyond the area sur veyed by these au thors. In or- 
der to complete the analysis of Mars ton et al.l ^1999), we study 
the molecular gas over a region ~ 86' x 86' around the star us- 
ing observations of the J=l— >0 transition of the CO and '^^CO 
molecules. Our aim is to investigate in more detail the morphol- 
ogy, kinematics and physical properties of the molecular envi- 
ronment of WR 16. To account for the properties of the ionized 
gas and dust we analyze available archival radio continuum and 
infrared (IR) data. The observations and databases used in this 
work are outlined in Sect. 2, the results are described in Sect. 3, 
and the discussion in Sect. 4. Concluding remarks are presented 
in Sect. 5. 



^ Radial velocities are referred to the local standard of rest (LSR) 



2. Observations and data bases 

2.1. CO data 

The CO J=l-^0 data were obtained using the 4-m NANTEN 
millimeter-wave telescope of Nagoya University. At the time 
we carried out the observations (April 2001), this telescope 
was installed at Las Campanas Observatory, Chile. The half- 
power beamwidth and the system temperature, including the at- 
mospheric contribution towards the zenith, were 2'. 6 (~1.7 pc 
at 2.3 kpc) and ~ 220 K (SSB) at 115 GHz, respectively. The 
data were gathered using the position switching mode with ab- 
solute reference positions; the points devoid of CO emission, 
used as reference, were available at the observatory and ob- 
served interspersed among the program positions. The spectrom- 
eter used was acousto-optical with 2048 channels with a fre- 
quency resolution of 20 kHz (equivalent to 0.055 km s"' at 1 15 
GHz) which provides a total coverage of ~1 13 km s"' . For in- 
tensity calibrations, a room-temperature chopper wheel was em- 
ployed (Penzias & Burrus 1973). An absolute intensity calibra- 
tion dUlich & Haaslll976l: iKutner & Ulichlll981 ) was achieved 
by observing Orion KL and p Oph East. The absolute radiation 
temperatures, T^, of Orion KL and p Oph East, which were both 
observed by the NANTEN radio telescope, were assumed to be 
65 K and 15 K, respectively (Moriguchi et al. 2001). A second- 
order degree polynomial was subtracted from the observations 
to account for instrumental baseline effects. The spectra were 
reduced using CLASS softwar^l The CO observations were 
centered at {l,b) = (281°04'47" -02°33'03") ((R.A., Dec)y20oo 
= 09''54"'53'' -57°43'38") with two areas covered with differ- 
ent samplings. An inner area of Al x Ab = 32'.4 x 32'.4 was 
sampled one beam apart, whilst a larger area of Al x Ab = 
86'. 4 X 86'. 4 was sampled every 5. '4 (two beamwidth). This 
technique allowed us to obtain two datacubes covering differ- 
ent areas around WR 16. The integration time per point was 16s 
resulting in a typical rms noise of ~ 0.3 K. 

2.2. CO data 

The radio astronomical group at Nagoya University has car- 
ried out a '^CO J=l— >0 survey of the Galactic Plane using the 
NANTEN telescope. In this work we use part of this unpublished 
survey that was kindly provided to us by the group. The beam 
size of the telescope at '^CO frequency (110 GHz) is 2. '7 and 
the typical rms noise is ~0.2 K. The frequency resolution is 40 
kHz (equivalent to 0.1 km s ' at 110 GHz). Observations were 
made in position switching mode and at a grid spacing of 4' 

2.3. Complementary data 

The CO and '^CO data are complemented with narrow-band Ha 
data retrieved from the SuperCOSMOS H-alpha Survey (SHS). 
The images have a sensitivity of 5 Rayleigh ( 1 Rayleigh = 2.41 X 
10"^ erg cm"^ s"' sr"'), and ~l" spatial resolution (Parke ret al.l 
2005). We also include radio continuum 4.85 GHz data retrieved 
from the Parkes-MIT-NRAO (PMN) Southern Radio Survey. 
The i mages have ~5 ' resolution and ~8 mJy beam"' total rms 
noise (ICondon et al.l[T99 3). To study the characteristics of the 
dust in the nebula we use infrared images at 60 pm and 100 
pm retriev ed from the high resolut ion IRAS images (HIRES) 
databas^ (iFowler & Aumannlll994 . 



3 htto://www.iram.fr/IRAMFR/PDB/class/class.html 
http://irsa.ipac.caltech.edu/applications/IRAS/IGA/ 
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3. Results and analysis of the observations 

3.1. Distribution and pliysical properties of tine CO and ^^CO 

We investigate the presence of molecular gas in the environs of 
WR 16 using several CO spectra obtained towards different po- 
sition around the star. After a careful inspection to the individual 
spectra, we note that the bulk of the CO emission is concentrated 

in the velocity range from 10 km s"' to — 3 km s"' . We use 

a mean radial velocity of 6 km s"' to make a rough estimate 

of the kinematical distance . Using the analytical fit to the circular 
galactic rotation model of iBrand & Bhtzl (11993.) along l^2%\° 
we estimate near and far kinematical distances of about 0.6 kpc 
and 2.8 kpc, respectively. Hence, the far kinematical distance is 
in agreement (within errors) with the adopted distance for WR 
16. This is indicative of a likely association between the molecu- 
lar gas detected in the mentioned velocity range and the WR star. 
We have to keep in mind, however, that non-circula r motions of 
up to 8 km s ' are know to exist in the Galaxy ( Burton & Gordon! 
[l978; Brand & Blitz 1993) which are of the same order of the 
velocity of the molecular feature. Although this might introduce 
large uncertainties in the kinematical distance determination it 
suffices for discerning between near and far kinematical dis- 
tance. 

We have analyzed the CO and '''CO data cubes in veloc- 
ity intervals of 0.2 km s ' whithin the mentioned velocity range. 
This strategy allowed us to identify three molecular features that 
are morphologically correlated with diff'erent regions of the ring 
nebula. The central velocities of these features are approximately 
at -5, -8.5, and -9.5 km s"' . With the aim to facilitate further 
analysis we shall refer to these features as component 1, com- 
ponent 2, and component 3, respectively. In the following, we 
describe and analyse the three features in detail. 

3.1 .1 . Component 1 

In Figs. 12^ and [2J3 we show the emission distribution of the 
'^CO and CO lines, respectively, in the velocity range from -7.2 
km s"' to -3.2 km s ' . Although the emission lines of CO and 
'■'CO trace the molecular gas of different densities, a molecular 
feature close to WR 16 stands out in both figures. This feature 
engulfs a low molecular emission region centered at the position 
of WR 16, surrounding the star from east to south. Fig.[2j3 also 
shows that the CO feature appears to be projected onto a more 
extended CO plateau having a relatively strong (~10 rms) emis- 
sion, that extends towards a large and bright molecular structure 
lying ~20' to the north and east, which is likely part of the galac- 
tic disk. 

In order to follow in more detail the morphology of compo- 
nent 1 and to compare with the Ha emission, we show in Fig. 
121; an overlay of the CO emission in the velocity range from 
-7.2 km s"' to -3.2 km s"' obtained with the 32'.4 x 32'.4 dat- 
acube, with the SHS emission of the ring nebula. Diff'erences in 
the contour lines between Figs.|2j5 and|2j; are due to the differ- 
ent sampling of b oth datacubes. This molecular feature was al- 
ready reported bv i Marston et al.l (119991) . referred to in that work 
as "molecular cocoon". 

We identify a number of clouds in the CO and '^CO emis- 
sions. They are labeled as al, a2, and a3 for the '^CO map, 
whilst for the CO image they are labeled as bl and b2. This 
identification is arbitrary and was made only to facilitate fur- 
ther analysis. Cloud al is very likely the '^CO counterpart of 
bl. On the other hand a2, and a3 are possibly the counterparts 
of b2 which appear as a single structure due to the worse CO 
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Fig. 2. Distribution of the mean CO emission integrated in the 
velocity range from -7.2 km s"' to -3.2 km s"' . The beam size 
of the CO observations is shown by a white circle in the lower 
left corner of each image. The position of WR 16 is marked by 
a plus sign at the center of each image, a) Emission distribution 
of the '^CO J=l— >0. The lowest temperature contour is 0.16 K 
(~ 5 rms). The contour spacing temperature is 0.08 K till 0.48 
K, and here onwards the contour spacing is 0.35 K. b) Emission 
distribution of the CO J=l— >0. The lowest temperature contour 
is 0.8 K (~10 rms) and the contour spacing temperature is 0.7 
K. c) Overlay of the mean CO emission in the velocity range 
from -7.2 km s"' to -3.2 km s"' (contours) and the Ha emis- 
sion (grey scale) in the central region. Contour levels are the 
same as in Fig 2b. 
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sampling. Assuming that the '^CO J=l— >0 line is in local ther- 
modynamic equilibrium (LTE), we follow the method described 
by Dickman ( 1978) to calculate some physical parameters such 
as the peak optical depth, t('^CO) and the molecular column 
density, A^('^CO). To use this method we assume that t('''CO) 
« 1, t(CO) » 1, and rexc(CO) = JexcC'^CO). Afterwards, we 
calculate the total hydrogen mass, M(H2), via 

M(H2) = (m,u„rVmH J] nN{n2)d^ (Mq) (1) 

where msun is the solar mass (~ 2 x lO-'^' g), fi is the mean 
molecular weight, which is assumed to be equal to 2.8 af- 
ter allowance of a relat ive helium abundance of 25% by mass 
dYamapuchietalJI 19991) . mH is the hydrogen atom mass (~ 1.67 
X 10"'^^ g), Q is the solid angle subtended by the CO feature in 
ster, d is the assumed distance expressed in cm, and A^(H2) is 
the H2 column density, obta ined using a "ca nonical" abundance 
N{U2)/NC^C0) = 5 X 10^ (lDickmanlll978l) . We also estimate 
the optical depth of the CO J=1^0 line from the '^CO J=l-^0 
line with t(CO)~(CO/'^CO)t(13CO), where CO/i^'CO is the iso- 
tope ratio (assumed to be ~62; iLanger & Penziaslll993l) . It is 
worth to point out that we adopted peak values of r*('^CO), 
which yields to upper limits for the line opacities, column den- 
sities, and molecular masses. 

The mass of the molecular gas can also be estimated us- 
ing the CO J=l— »0 line via two methods: (a) Using the relation 
between integrated emission of the CO J=l— >0 line, /(CO) - 
J r^(CO) dv, and molecular hydrogen column density, A^(H2), 
and (b) the virial theorem. For the first method we use 

A^(H2) = X X /(CO) (cm"2) (2) 

The X-factor Ues in the range (1-3) x 10^° cm"^ (K 
km )~\ as estimated by the virial theorem and y-ray emission 
(iBloemen et al. 1986; Solomon etal. 1987: Be rtsch et al .1119931: 
[Pigel et a l. 1996; Strong & Mattox 1 996). In this paper we adopt 
X = 1.6 X 10^" cm-2 (K km s ' )-' dHunter et al.llI997h . For the 
second method, considering only gravitational and internal pres- 
sure and assuming a spherically symmetric cloud with a r"' den- 
sity distribution, the molecular mass can be determined from 

M = 190/J(Avcid)^ (Me) (3) 

(iMacLaren et al.l[l988h . where R- V^cioud/^ is the effective ra- 
dius in parsecs, and Avdd is defined as the FWHM line width of 
the composite profile derived by using a single Gaussian fitting. 
The composite profile is obtained by averaging all the spectra 
within the area of the cloud (Adoud)- The uncertainties involved 
are about 40 %, mostly due to the estimates of R. In Table[T]we 
summarize some of the physical parameters detailed above, ob- 
tained for the molecular clouds constituting component 1 . The 
values of A^(H2)x and M(H2)x are the hydrogen column den- 
sity and molecular mass, respectively, determined using Eq|2] 
whilst A^('^CO)lte and M(H2)lte are the '^CO column density 
and molecular mass obtained assuming LTE. The virialized mass 
is given by Mvii . In Table [T] we also list the thermal line widths 
for bl and b2, Av±, obtained assuming Tkin = Texc (LTE). To 
derive the total molecular mass of component 1 from the '^'CO 
data we added up the values obtained for al, a2, and a3, which 
yields to a total mass M(H2)lte ~ 900 Mq. Using the same tech- 
nique with the values obtained for bl and b2 from the CO data, 
we derive total molecular masses M(H2)x ~ 5500 and Myir 
^ 16000 Mo. 

It is worth to point out that the physical parameters quoted in 
Table [T] are only derived for component 1, since is the only one 
entirely detected at '^CO emission and morphologically defined 
without ambiguity. 



3.1.2. Component 2 

In Fig. [3] we show the spatial distribution of the CO emission in 
the velocity interval from -9.5 km s"' to -7.2 km s"' , overlaid 
on the Ha emission of the nebula. The molecular gas in this 
velocity interval shows a morphological correspondence with 
the outer Ha ring, especially at the intense optical emission at 
(R.A., De c.)72ono ^ (9^56'"15V 57°51'). This feature was not de- 
tected by iMarston et al.l ( 1199 9^ due to the limited area covered 
by their CO data. It can be noted that this feature extends well 
beyond the area depicted in Fig. |3] which suggests that is part of 
a larger structure, probably the previously reported extended CO 
plateau lying aside component 1 (see Fig.|2j)). 




9 57 00 56 00 55 00 54 00 53 00 



Right Ascension (J2000) 

Fig. 3. CO integrated intensity in the velocity range from -9.5 
km s"' to -7.2 km s"' (contours) superimposed onto the Ha 
emission (grey scale). The contour levels start at 0.8 K and the 
contour spacing is 0.55 K. 

In Fig.|4]we show the the CO emission distribution of com- 
ponent 2 (grayscale) superimposed on that of component 1 (con- 
tours). A spatial separation between both components can be ob- 
served, with component 2 adjacent to the eastern border of the 
outer Ha ring, and component 1 placed almost at the interface 
between the latter and the Ha emission. The kinematic of the 
molecular gas was analyzed by position-velocity maps across 
selected strips. In the right panels of Fig. |4] we show the position- 
velocity diagrams in the directions sketched by the dotted lines 
in the left panel of the figure. Slices a to c show that the peak 
emission of both features are discernible and spatially separated 
by about ~ 7'. Both slices also show the presence of a weak 
bridge of CO emission which connects the western border of 
component 1 with the eastern side of component 2. Although 
this bridge might suggest a physical link between both features, 
no hints of a velocity gradient is detected. On the other hand, 
slice c shows a stronger CO emission bridge between the two 
features. However, the peak emission of both features are well 
separated by ~ 6'. Probably, the selected direction of slice b, 
along the stronger region of component 2, and the weaker region 
of component 1, makes their emissions to look merged, although 
this may be due to a projection effect. The morphology and kine- 
matical behavior of component 2 suggests that the presence of 
this feature is not consequence of a mere velocity gradient, but 
to the presence of a second molecular structure. As suggested 
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Fig. 4. Left panel: Comparison between component 1 (contours) and component 2 (grey scale). The location of the slices, labeled as 
a, b, and c, are sketched by the dotted lines. The cross at the center indicates the position of WR 16. Right panels: Position-velocity 
maps obtained for strips a, b, and c. Contour levels start at 1.9 K and the contour spacing is 1.4 K. The loci of both components are 
indicated by the black arrows. 



before, component 2 may be a small part of a larger molecular 
structure (plateau), whose extension likely exceeds the area de- 
picted in Figs. [3] and |4] The double-layer morphology reported 

above for the molecular gas in the velocity range from 10 

km s"' to — 3 km s ' could be associated to some shock fronts 
and mass-loss rate episodes which have been acting upon the 
surrounding ISM, as p redicted by the hydrodynamical models 
of Garci a-Segura et al.l (1996a b). Characteristics times of this 
features are further discussed in Sect. 4.2. 

3.1.3. Component 3 

In Fig.|5]we show an image ~ 16' x 16' in size centered at the 
position of WR16 that shows the CO emission in the velocity 
range from -9.59 km s ' to -9.48 km s ' (two channel maps), su- 
perimposed onto the Via emission of the inner ring nebula. 

Even though the intense CO clump located at 
(R.A., Dec.)/20oo = (9''55"'30'', -57°47'30") stands out, a 
fainter ring-shaped molecular emission (indicated in the fig- 
ure as component 3) surrounding the position of WR 16 is 
noticeable. The star is seen projected onto a well-defined 
minimum of the CO emission whilst the molecular gas ~ 3' 
north of the star is coincident with the brightest Ha emission 
of the inner ring. The good morphological correspondence 
between the CO and the Ha emission may be indicative of 
a mutual physical association between the molecular and the 
ionized gas . It is worth t o mention that this feature was not 
detected bv lMarston et alj {1999), probably due to sensitivity or 
velocity resolution. However, we are cautious about reporting 
component 3 since the area is small and the CO observations 
are not fully-sampled. Furthermore, we note the short velocity 
interval in which component 3 is noticeable. Further high-spatial 
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Fig. 5. Distribution of the mean CO emission in the velocity 
interval from -9.59 km s"' to -9.48 km s"' (contours), superim- 
posed on the Ha emission (grey scale). The beam size of the 
CO observations is shown by a white circle in the lower left cor- 
ner The lowest temperature contour is 0.5 K (~ 4.5 rms). The 
contour spacing is ~0.65 K. 



resolution observations with instruments like APEX may help 
to clarify the nature and characteristics of this feature. 
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Table 1. Main physical parameters of the identified molecular 
clouds 



' 'CO clouds 


Parameters 


al 


a2 


a3 


Center position: 

K-A(J2000) I, ) 


09 56 30 


09 54 30 


09 53 00 


DeC(j2000) (° "') 


-57 37 30 


-57 59 30 


-57 57 00 


r-R(Peak) (K) 


0.45 


0.53 


0.32 


tC^CO) 


-0.09 


-0.11 


-0.07 


T(CO) 


-1.7 


-2.5 


-1.1 


£1(10^' ster) 


7.4 


9.7 


5.1 


A'(''CO)lte(10'= cm-2) 


<0.7 


< 1.1 


<0.15 


MHz)™ (10^" cm-2) 


<3.4 


<5.3 


<0.8 


M(H2)lte (10^ MO) 


- 2.8 


-5.6 


-0.5 


CO clouds 


Parameters 


bl 


b2 




Center position: 








KA(J2000) (, ) 


09 56 00 


09 54 30 




DeC(j2ooo) (° ' ") 


-57 43 00 


-58 02 00 




r*R(Peak) (K) 


3.9 


4.1 




T-exc (K) 


7.1 


7.3 




r*R(average) (K) 


2.1 


2.2 




/(CO) (Kkms-') 


- 6.3 


- 6.6 




£1(10"'' ster) 


22.8 


24.3 




A'(H2)x (10^" cm-2) 


10.1 


10.6 




M(H2)x (10^ MO) 


2.6±1.1 


2.9±1.3 




Avcid (km s"' ) 


3.3 


4.1 




Av,i, (km s"' ) 


0.10 


0.11 




Ripe) 


2.9 


3.1 




(10^ MO) 


-6.1 


-9.9 





3.2. Radio continuum emission 

To show the distribution of the ionized gas, we display in Fig. 
|6]the 4.85 GHz radio continuum emission image obtained from 
the PMN survey, overlaid onto the SHS Ha emission image 
of the nebula. As expected, both images show a relatively good 
morphological correspondence, since both emissions are indica- 
tive of the existence of ionized gas. Four radio continuum fea- 
tures are identified within the area observed at CO frequencies. 
Feature 1 peaks at (R.A., Dec.)y20()0 ~ (9''54'"40'',-57°40'56") 
and has an arc-like morphology almost coincident with the 
brightest section of the inner ring nebula. Feature 2, peaking 
at (R.A., Dec.)y20oo ~ (9''56'"40',-57°51'30"), is placed along- 
side a relatively strong optical structure. The peak emission of 
Feature 3 (R.A., Dec.)y20oo ~ (9''55'"00',-57°35'00") is coinci- 
dent with the brightest section of the outer ring nebula. Unlike 
the other radio continuum features. Feature 4 is seen projected 
onto a weak optical emission. Its location and morphology sug- 
gests that this feature might be the radio continuum counterpart 
of the south-western section of the outer ring nebula. 

In Table |2] we summarize some parameters of the ionized 
gas estimated from the radio continuum flux densities at 4.85 
GHz (54.85). Adopting an electron temperature of 10"^ K, assum- 
ing co nstant electron densities, and using the spherical models 
of iMezger & HendersonI d 19671) for features 2 and 3, and their 
cylindrical model for features 1 and 4, we obtained the rms elec- 
tron density (we) and the ionized mass (Mion) for each feature. 
For 5 4.85 we estimate an eiTor of about 20%. Since no clumping 
factor was taken into account, the obtained values of «e and Mion 
result in lower and upper limits, respectively. To calculate the 
number of Lyman continuum photons needed to keep the cur- 



in the environs of the star WR 16 
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Fig. 6. Overlay of the radio continuum image at 4.85 GHz (con- 
tour lines) and the Ha image (grey scale). The lowest contour 
line is 28 mJy beam"' (~ 3.5 rms), and the contour spacing is 50 
mJy beam 

rent level of ionizat i on in the ring nebulae (A^Lyc), we use the 
relation o f IChaissonI (119761) . 

A^Lyc = 0.76 X 10^^ T^°-^^ v" ' Sy cp-, (4) 

where T4 is the electron temperature in units of lO'* K, v is the 
frequency in units of GHz, d is the distance in kpc, and 5v is 
the total flux density in Jy. Summing up the values obtained for 
all the radio continuum features we obtain A^Lyc~l x s"'. 
Considering that the number of Lyman continuum photons emit- 
ted by a WN8 star is - L25 x lO^*^ s"' (Crowther 2007), we 
therefore conclude that WR 16 is by far capable of maintaining 
the observed ionization level. 

Table 2. Main parameters of the ionized gas in features 1, 2, 3, 
and 4 



features 


1 


2 


3 


4 


54.85 (mJy) 


-48 


-24 


-95 


-50 


He (cm"') 


16 


10 


7 


8 


Mion (Mo) 


4 


3 


20 


9 


A'Lyc (lO-^^seg-') 


0.2 


0.1 


0.5 


0.2 



4. Discussion 

4. 1 . The moiecular mass of component 1 

In Sect. 3.1 we have derived the mass of the molecular clouds 
constituting component 1 making use of the CO and '^CO 
data by three diff'erent methods. In spite of their uncertainties, 
all three methods provide mass estimates for this feature that 
strongly suggest that mostly of the molecular gas has a non- 
stellar origin. We note that the total LTE mass obtained in this 
work (-1200 Mq) is substantially larger than that obtained by 
Ma rston et al. l (IT999l) . although this difference may be due to the 
limited area covered by these authors. 



6 



N. U. Duronea et al.: Carbon monoxide in the environs of the star WR 16 



According to the obtained values a remarkable result is no- 
ticed when comparing all three methods, namely Mvir > M(H2)x 
> M(H2)lte- Regarding to the total virialized mass, the derived 
value (~ 16000 M©) was obtained assuming that only the grav- 
itational and internal pressure support the stability of the cloud 
ignoring the effect of external pressure. Mass estimation meth- 
ods such as the Virial and CO integrated emission masses depend 
on the line width, which are highly sensitive to cloud kinematics. 
Hence, comparing masses estimated by these methods enables 
some analysis of turbulence and kinematics of the molecular gas. 
From Eqs.|2]and[3]it can be noted that MyJM(H2)x °^ Avcid/^- 
Therefore, molecular clouds having high line widths should ex- 
hibit large virialized mass in comparison with its integrated 
mass. This could be the case of component 1, since bl and b2 
exhibit line widths significantly higher than the predicted ther- 
mal width (see Table [U. This is not odd since the environs of a 
massive star is highly disturbed by strong stellar winds and vio- 
lent shock fronts. The high value of Mvii with respect to M(H2)x 
indicates that component 1 has too much kinetic energy and is 
not in virial equilibrium. Very likely, the molecular gas of com- 
ponent 1 is supported by the external pressure provided by the 
warm and low density gas that is ionized by the strong UV ra- 
diation field of the WR star. The pressure needed to sustain the 
molecular gas could be provided by the ionized gas of the ring 
nebula trapping component 1 . 

The discrepancy between M(H2)x and M(H2)lte is harder 
to explain. The X-factor was obtained empirically towards sev- 
eral molecular clouds in the Galactic disk, and is believed to 
be accurate to a factor of 2, when averaged over a large num- 
ber of clouds. However, this factor can vary under conditions 
differing from those of typical Galactic clouds, namely density, 
age, turbulence, cosmic ionization rate, incident radiation field, 
and metalHcity (iBell et al.ll2006l:lPineda et"ani2008l;IShettv et alj 
1201 1 ). The last two factors may have played a major role in the 
environs of the star In fact, these two factors are intimately re- 
lated since the formation of CO is sensitive to both, the amount 
of available C and O, and the UV radiation field responsible for 
photodissociating molecules and regulating the CO abundances. 
Increasing metallicity could strongly decrease the X-factor be- 
low the canonical value, whilst a stronger radiation field could 
slightly raise it ( Bell et al.l 12006). The chemical abundance of 
the inner ring nebula around WR 16 suggests that the close sur- 
roundings of the WR star have been processed and enriched by 
molecular ejecta orig inated during previous LBV or RSG phases 
(iMarston et al.lll999l) . which very likely dominates the produc- 
tion of new CO despite its depletion due the ionizing radiation 
field of the star. Then, component 1 may be actually composed 
by a mix of stellar ejecta and interstellar molecular gas. This 
could mean that the X-factor used in Eq. [2] (1-6 x 10^" cm"^ 
(K km s"' ) ') is high for this case, which originates some dif- 
ference between M(H2)x and M(H2)lte- In fact, this might also 
explain the high ratio r*R(CO)/r*R('^CO) observed for com- 
ponent 2 (~ 9; see Table [TJ, which is almost twice the average 
Galactic value. It is worth to point out that turbulent motions pro- 
duce a slig ht drop in the va lue of X, although this effect is not 
significant (iBell et al ]l2006h . Nevertheless, we keep in mind that 
'^CO LTE-based mass estimates are also subject to their own 
uncertainties and may underestim ate true column density by a 
factor up to 7 (iPadoan et alj|2000h . 

The study of higher rotational levels of the CO could shed 
some light on the physical and excitation conditions of the 
molecular gas in the region. Further, non-LTE analysis of these 
lines will help to obtain more robust estimates of some physical 



parameters, such as the column density, mass, kinetic tempera- 
ture, and opacities. 



4.2. Dynamics of the molecular gas 

In previous sections we have detailed the characteristics of the 
molecular gas components projected onto the multiple ring neb- 
ula around WR 16. The good morphological correspondence be- 
tween the CO and Ha emissions is suggestive of a physical as- 
sociation between the molecular and the ionized gas. 

The morphology of the molecular gas in the velocity range 
from — 10 km s"' to -3 km s"' suggests that the ambient molec- 
ular gas in the vicinity of WR 16 could have been swept up by 
the action of shock-fronts and the strong stellar winds, giving 
rise to the features shown in Fig. |4l Assuming that the winds 
of the star are spherically symmetric, and that component 1 and 
component 2 are two-dimensional projections of spherical ex- 
panding shells, a connected structure with the velocity should 
be manifested in the datacube, first as a blueshifted pole, then 
into a growing and decreasing circular ring, and finally to a red- 
shifted pole, with the powering star placed at the center or close 
to it. However, according to our data none of these kinemati- 
cal features are noticed, since no hints of a velocity gradient is 
detected. As we have pointed out, component 2 is very likely 
part of a larger molecular emission plateau that extends to the 
east (see Fig. |2]i. We infer that WR 16 is placed outside of, 
or in an eccentric position with respect to this plateau, and is 
blowing the molecular gas from a side (face-on). Evidences of 
this kind of molecular structures have been found in the Galaxy 
(see lArce et al.ll20Tll) . The location of the brightest region of 
the outer ring nebula, almost in the opposite direction of com- 
ponent 1 and component 2, suggests that this ring might be the 
remains of the parental molecular gas after been ionized by the 
WR star and/or its progenitor The molecular clump located at 
(R.A., Dec.)y2000 ~ (9''54",-57°33') could be a remnant of the 
parental molecular gas (see Fig.[3]i. 

As mentioned in Sect. 3.1, the stratified location and kine- 
matics of component 1 and component 2 (and probably compo- 
nent 3) is highly correlated with the multiple layer morphology 
of the nebula, and may be indicative of an evolution of the ISM 
associated with successive shock-fronts, each one related with 
different evolutive episodes of the predecessors of WR 16. The 
different velocities of each component could have been produced 
by violent events related to the star evolution and the probable 
occurrence of mass eruptions, which combined with the strong 
stellar UV radiation field, have produced the Ho- morphology 
and the observed kinematical features. Component 2 seems to 
be the more extended. Its border is adjacent to the eastern border 
of the outer Ha ring. Probably this is the oldest one and could 
represent the main-sequence signature of the star progenitor (see 
below). Component 1 is placed almost at the interface between 
component 2 and the Ha emission. It could be related to an in- 
terme diate phase prior to the actual WR stage, such as RSG or 
LBV dMarston et al.lll999h . Finally, the location of component 3 
is internal to the other two components, and its origin could be 
possible associated with a youngest event. 

The mentioned scenario is rem arkably similar to that de- 
picted bv lRizzo et al ] (l2001bll200 ?) forNGC 2359. In that case, 
however, the main "fingerprint" of the previous 0-stage was re- 
vealed by the presence of a Hi shell expanding at 12 km s ' that 
shocked the molecular gas accelerating it to a velocity of ~ 12-13 
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km s . In order to test the proposed scenario a rough estimate 
of the dynamical age (fdyn) of component 2 was made using 

fdyn = 0.5 X 10" /^/Vexp (yr), (5) 

dHowarth & LamersI [19991) . Here, R is the averaged radius of 
component 2 (~ 9 pc at 2.3 kpc), and Vgxp is the expansion veloc- 
ity. For Vexp we consid er a conservative range of values from ~5 
km s"^ dPuronea et al.ll2012h to -12 km s^^ dRizzo et al.ll2003h . 
Then, the obtained values of fdyn are in the range 4 - 9 x 10^ 
yr. If we consider t he 35 M0 model of stellar evolution of 
iGarcia-Segura et all dl996bi) (O^RSG^WR), the duration of 
the O phase (4.5 X lO" yr) is larger than fdyn by a factor of ~ 5 - 
10. On the other ha nd, taking into account the 60 Mo model of 
iLanger et alJ(fT994l ) (O^pre-LBV^LBV^WR) the duration of 
the O phase is larger only by a factor ~ 2 - 4. In any case, these 
differences could be explained in terms of the small solid angle 
subtended by component 2 around WR 16 (see Fig.O, which is 
not able to be impelled by the totality of the stellar wind power 
This scenario allows one to deal with the hypothesis by which 
the parental molecular cloud around an evolved massive star is 
destroyed and/or evacuated after star formation. It is worth to 
point out that we have searched for Hi shells around WR 16, 
using available surveys, having negative results. This reinforces 
the proposed scenario of component 2 being blown-up by the 
O-star progenitor 

The morphology and location of component 1 might lead one 
to assume that this feature is composed of molecular gas ejected 
by the star in a previous evolutionary phase. However, the molec- 
ular mass estimated in Sect. 3.1 led us to conclude that this com- 
ponent might be rather composed by a mix of shocked inter- 
stellar molecular gas and stellar ejecta from intermediate stages 
(RSG or LBV), proba bly as a result of th e encounter of two 
winds (as suggested bv lGarcia-Segura et alJ[T99 6a b) impinging 
onto the interstellar molecular gas. The most likely scenario is 
the WR wind encountering an already stopped RSG/LBV wind 
in contact with the molecular cloud, and being abruptly deceler- 
ated giving rise to the observed shocked gas. The stellar ejecta 
has been probably accum ulating at th e inner edge of compo- 
nent 1, as proposed by IMarston et alJ (11999'), where the WR 
wind is impacting and the UV field is ionizing, giving rise to 
the inner Ho' ring. The velocity difference among component 
1 and component 2 (2 - 3 km s"' ) would give a rough esti- 
mate of the velocity of the shock, although this is very likely 
a lower limit since we are observing projected radial velocities 
from a face-on geometry. The location of component 1 entails 
that the wind of the previous O-stage have not fully interacted 
with the ISM and could have partially run through the molec- 
ular gas. Indeed, this might be possible owing to the relatively 
low densi ty of the O-stage w inds. This argument was also dis- 
cussed bv lRizzo et alj d2003i) . and this scenario implies that only 
a fraction of the mechanical energy injected by the winds of the 
previous O-stage is converted into kinetic energy of the expand- 
ing molecular gas, which is in line with theoretical predictions 
dKoo & McKedll992l:lAr5iu3l2007h . 

4.3. The inner ring 

In this section we will briefly discuss so me aspects of the in- 
ner region of the nebula. In their Fig. 8. IMarston eTal] d 19991) 
show the continuum 60 //m emission distribution around the star 
WR 16. The IRAS map reveals a bright shell that has an ex- 
cellent spatial correlation with the inner ring nebula (no other 
external IR features were reported in the paper). This shell is 



al so detected in the HI RES IRAS 25 yum emission (see Fig. 9 
of Marst on et alJ[i999l) . We have analyzed the dust continuum 
emission at 60 /urn and 100 fim from IRAS high resolution im- 
ages (HIRES) to study the warm du st counterparts of the inner 
ring nebula. Following Draind d 19901) the dust temperature of the 
shell was obtained using 

Td = 95.94/ln[l.667P+">(Fioo/i^6o)] (K) (6) 

where Fioo and Fgo are the 100 fim and 60 fim IRAS fluxes. 
The parameter n (assumed to be 1.5) is related to the absorption 
efficiency of the dust. We obtained a dust temperature of about 
Td ~ 28 K, which is re markably simil ar to that found in the ring 
nebula NGC 2359 (Riz zo et alji2001ah . We also derived flie dust 
mass of the shell using 

= 4.81 X 10"'2 ^100 'i^(e"*^ '*^^^''-l) (Mq) (7) 

(Tokunaga 2000"), where d the distance in pc, and Fioo is given 
Jy. We obtain a dust mass of ~ 0. 1 Mowith an uncertainty of 
the order of ~ 45 %. By assuming a standard gas-to-dust ratio 
of ~ 100, and taking into account an ionized gas of ~ 4 M0(see 
Sect. 3.2), a molecular gas component of ~ 6 M© in mass is 
predicted to be associated with the dust shell around the WR 
star, which could explain the presence of component 3. 

In Fig.|7]we show an overlay of the mean CO line emission 
in a velocity range -9.59 km s"' to -9.48 km s"' (component 
3) onto the HIRES 60 fim emission which has been smoothed 
down to a common beam of 2. '7. A look at this figure shows 
that the 60 /jm and molecular emission are asymmetrical, 
having most of their emission toward the north western section. 
The same behavior is observed for the optical and the radio 




9 55 45 30 15 55 00 45 30 15 00 



Righ Ascension J(2000) 

Fig. 7. Overlay of the mean CO fine (contours) in the velocity 
range from -9.59 km s ' to -9.48 km s"' , and HIRES 60 yum 
emission (grey scale) smoothed down to a beam of 2'. 7. Dashed 
and soli d arrows indi cate the direction of the proper motion de- 
rived bv lMoffat et al.l ([l998|,il999) and in this work, respectively. 

continuum emissions (see Figs.[T]and|6]l. Using pro per rnotions 
obtained from the Hipparcos astrometry satellite (! ESAl[T997h 
Moffat etal. (1998, 1999) suggest that WR 16 is a candidate 
star to be a runaway. The direction of the peculiar tangentia l 
proper motion of WR 16 (according to lMoffat et al.lll998lll999l) 



8 



N. U. Duronea et al.: Carbon monoxide in the environs of the star WR 16 



is shown in Fig. [T] by a dashed arrow. Follo wing Moffat et"an 
(Il998lll999l) but using the Tycho-2 Catalogue (Hos et alJl2000h 
we derived for WR 16 a new direction of the peculiar tangential 
velocity (fii^pec, f^b,pec) with respect to its local ISM. Adopting 
observed proper motions of //„cos((J) = (-6.8 + 1.4) mas yr"' 
and = (5 ± 1.4) mas yr"' we derived (tii,pec, f^b,pec) = 
(-12.1 + 1.1) mas yr"', (-3.5 + 2.6) mas yr"'. The obtained 
direction of the peculiar proper motion computed in this way 
is depict ed in Fig. U\ by a solid arrow. Our results confirm 
those of iMoffat et alT(ll999 ') and very likely suggest that the 
strengthening in the IR, radio-continuum, optical, and probably 
molecular emission may be caused by the motion of the WR 
star with respect to the ISM ("bow-shock"). 

From the explained above, it becomes clear the need of better 
evolutionary models which may explain the existence of dust and 
molecular gas at such close distance to a WR star. However, as 
pointed out in Sect. 3 . 1 , it is important to stress that the existence 
of component 3 is still to be confirmed by high-resolution and 
sensitive observations. 



5. Concluding remarks 

We have performed a study of the molecular gas in the mul- 
tiple ring nebula associated with the WN8 star WR 16. In 
this work, NANTEN CO J=1^0 and ^^CO J=1^0 observations 
were used, together with 4.85 GHz radio continuum data re- 
trieved from the PMNRAO Southern Radio Survey, and IRAS 
(HIRES) archival data. 

The new CO and '^CO data allowed us to map the whole 
molecular gas surrounding the star WR 16. We reported two 
molecular features at velocities of -5 km s"' (component 1) and 
-8.5 km s"' (component2). For component 1 we derive a molec- 
ular mass > 900 Mq which indicates that is mostly composed by 
interstellar gas. We also find that Mvk > M(H2)x > M(H2)lte 
which suggests that this feature is probably binded by the exter- 
nal pressure of the ionized gas, and is enriched by stellar molec- 
ular ejecta. The morphological and kinematical stratification of 

the molecular gas in the velocity range from 10 km s~' to -3 

km s ' is indicative of an evolution of the surrounding molecular 
gas with successive events associated with the evolution of the 
WR star. Based on our findings, it is suggested that WR 16 could 
have been located close to the edge, or outside of a large molec- 
ular structure. The winds of the progenitor O-star could be re- 
sponsible of the formation of component 2, ploughing part of the 
molecular gas surrounding the star The rest of the molecular gas 
was first impacted by the winds of previous RSG or LBV phase, 
and later by the WR winds, which originated a shock front, and 
the kinematical features observed in the CO emission. 

The 4.85 GHz emission shows a relatively good morphologi- 
cal correspondence with the optical nebula. The derived electron 
densities are about 16 cm"-* for the inner nebula, and 8-10 cm""* 
for the outer nebula. The total ionized gas is ~36 Mp. 

Th e dust shell around WR 16 reported by iMarston et alj 
(11999^ was studied using HIRES IRAS 60 jum and 100 yum. We 
determined a dust temperature and mass of about 28 K and 0. 1 
Mq, respectively. We also reported a third molecular feature at 
-9.5 km s"' (component 3), having a good morphological corre- 
lation with the the dust continuum emission at 60 jum and with 
the Ha emission, which has to be confirmed with high resolution 
observations. The asymmetry of the discovered structures may 
also be qualitatively explained by the motion of the massive star 
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